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Abstract 

We generalize Hagedorn's statistical theory of momentum spectra 
of particles produced in high-energy collisions using Tsallis' formalism 
of non-extensive statistical mechanics. Suitable non-extensive grand 
canonical partition functions are introduced for both fermions and 
bosons. Average occupation numbers and moments of transverse mo- 
menta are evaluated in an analytic way. We analyse the energy depen- 
dence of the non-extensitivity parameter q as well as the g-dependence 
of the Hagedorn temperature. We also take into account the multiplic- 
ity. As a final result we obtain formulas for differential cross sections 
that are in very good agreement with e + e~ annihilation experiments. 
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1 Introduction 



35 years ago Hagedorn has written a paper of fundamental importance Jl| 
which nowadays is very often cited in various contexts of statistical me- 
chanics, particle physics, and string theory. In this paper he developed a 
statistical description of momentum spectra of particles produced in collider 
experiments. At a fundamental level, of course, the underlying theory for 
this is quantum chromodynamics. However, for the hadronization cascade 
phenomenological models have to be used in addition to QCD. The complex- 
ity inherent in this process is so immense that usually one has to perform 
Monte Carlo simulations in order to explain the experimental measurements 
of cross sections. 

Hagedorn's successful approach was to see this problem from a thermo- 
dynamic point of view (of course QCD was not known at the time he wrote 
the paper). He devoped a theory that is nowadays regularly in use do de- 
scribe e.g. heavy ion collisions |2|, |^, M. The basic assumption is that in the 
scattering region the density of states grows so rapidly that the temperature 
cannot exceed a certain maximum temperature, the Hagedorn temperature 
To. This is similar to a first oder phase transition of e.g. boiling water, just 
that the Hagedorn temperature describes 'boiling' nuclear matter at about 
T ~ 200 MeV. This state of matter is closely related to what is nowadays 
called quark gluon plasma. For some theoretical approaches to transverse 
momentum spectra, see e.g. || ||, 0. As already mentioned, the Hagedorn 
phase transition is also of fundamental interest in string theories || || |10| . 

While for collider experiments (e.g. e + e~ annihilation) the Hagedorn 
theory yields a good description for relatively small center of mass enegies 
(E < 10 GeV), it fails at large energies. Indeed, Hagedorn's approach pre- 
dicts an exponential decay of differential cross sections for large transverse 
momenta, whereas in experiments one observes non-exponential behaviour 
for large energies (E > 10 GeV) |JTT], [12], [13| . For this reason one usually re- 
stricts the comparison of the experimental results at higher energies to Monte 
Carlo simulations, which reasonably well reproduce the experimental data. 
But apparently, for a deeper understanding from a statistical mechanics point 
of view there is the need to generalize Hagedorn's original ideas. 

In this paper we show that it is possible to extend Hagedorn's theory in 
such a way that it correctly describes the experimental findings at large en- 
ergies as well. The basic input for this is the recently developed formalism of 



non-extensive statistical mechanics [|rjL [TR (TBI, a generalization of ordinary 
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statistical mechanics suitable for multifractal and self-similar systems with 
long-range interactions. The formalism has been introduced by Tsallis 12 



years ago ]TJJ and has in the mean time been shown to be very successful to 



describe e.g. systems exhibiting turbulent behaviour [17, 18, [L9[ or HamiL 



tonian systems with long-range interactions ||20|1 . Several other physical ap- 
plications were described in — [p5|| . Very recently, evidence was provided 
that the new formalism has applications in particle physics as well. Wal- 
ton and Rafelski [[26j studied a Fokker Planck equation describing charmed 
quarks in a thermal quark-gluon plasma and showed that Tsallis statistics is 
relevant. Wilk and Wlodarczyk [^7j provided evidence that the distribution 
of depths of vertices of hadrons originating from cosmic ray cascades, as mea- 
sured in emulsion chamber experiments, follows Tsallis statistics. Alberico, 
Lavagno and Quarati |2£| successfully analysed heavy ion collisions using 
Tsallis statistics. Bediaga, Curado and Miranda showed that differential 
cross sections in e + e~ annihilation experiments can be very very well fitted 
using Tsallis statistics. In the present paper we will work out the theory 
underlying this. 

To effectively describe the complex QCD interactions and hadronization 
cascade by a a thermodynamic model, it seems reasonable to use the Tsallis 
formalism (which contains ordinary statistical mechanics as a special case), 
since this formalism is specially designed to include selfsimilar systems and 
systems with long-range interactions. What happens in a collider experiment 
is expected to fall into this category, since QCD forces are strong at large 
distances. Moreover, the self-similarity of the scattering process was already 
recognized by Hagedorn, who described the various possible particle states 
as a 'fireball' and who defined a fireball as follows f50fl: A fireball is 



a statistical equilibrium of an undetermined number of all kinds of 
fireballs, each of which in turn is considered to be... 
(back to *) 

Clearly, nowadays we would call this a self-similarity assumption. 

In this paper we will combine Hagedorn's and Tsallis' approach. We will 
introduce a suitable generalized grand canonical partition function. The 
predictions following from this generalized thermodynamic model can be 
evaluated analytically (at least in a certain approximation). We will also 
take into account the multiplicity and suggest a concrete dependence of the 
non-extensitivity parameter q on the energy E. The result will be a con- 
crete formula for differential cross sections of transverse momenta pt where 
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no free fitting parameters are left. Our formula turns out to be in excel- 
lent agreement with experimental measurements in e + e~ annihilation exper- 
iments. This turns out to be true for the entire range of center of mass 
energies that have been probed in experiments. Our analytical formulas are 
also in good agreement with the curves obtained from Monte Carlo simula- 
tions, thus suggesting that the complexity inherent in these simulations is 
effectively reproduced by the simple thermodynamic model considered here. 

This paper is organized as follows. In section 2 we shortly review the 
main results of Hagedorn's theory. In section 3 we shortly review Tsallis' 
theory. In section 4 we combine both. We will write down the generalized 
grand canonical partition function, including both fermions and bosons. In 
section 5 we will consider a large pt approximation, which allows for a simple 
analytical treatment. We will numerically show that this approximation is a 
good one even for relatively small values of the transverse momentum. We 
derive the relevant form of the probability density and calculate all moments 
of transverse momenta. In section 6 we investigate the (weak) dependence 
of the Hagedorn temperature T on the non-extensitivity parameter q. In 
section 7 we analyse the energy dependence of the parameter q. In section 
8 we take into account the multiplicity and derive the final formula for the 
differential cross section. Finally, in section 9 we compare our theoretical 
results with the experimentally measured results of the TASSO and DELPHI 
collaboration. Our concluding remarks are given in section 10. 



2 Hagedorn's theory 

Hagedorn's theory IJ essentially predicts that the differential cross section 
in a scattering experiment is given by 



1 da 



CPT 



DC 



dx e -/V* 2 +^ 2 (1) 



<T dp T JO 

Here pr is the transverse momentum, /i = \jp\ + m? is the transverse 
mass, the integration variable x stands for the longitudinal momentum, 
(3 = l/(kT ) denotes the inverse Hagedorn temperature, and c is some con- 
stant. The factor e-P\/ x2+ v 2 is immediately recognized as a Boltzmann factor 
with energy given by the relativistic energy-momentum relation. In the fol- 
lowing, we will set the Boltzmann constant k equal to 1. 
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The Hagedorn temperature T is independent of the center of mass energy 
E of the beam and typically has a a value of about 100-300 MeV. The physical 
idea underlying Hagedorn's approach is that To describes a kind of 'boiling 
temperature' of nuclear matter, which cannot be further increased by external 
energy transfer. Any further increase of energy just produces new states of 
particles rather than an increase of temperature. 

The integral in eq. (P can be further evaluated to yield 

where K\ is the modified Hankel function. For p? large compared to both 
To and m one obtains from the asymptotic behaviour of the Hankel function 
the approximate formula 

I^pfe"^. (3) 
a dp T 

One sees that the differential cross section decays exponentially for large 
values of pr- While this exponential decay is indeed observed for collider 
experiments with relatively small center of mass energies (E < 10 GeV), 
clear deviations from exponential decay have been observed at higher en- 
ergies |TT|, 0, |I3[ . Here the Hagedorn theory is not valid any more and 
the dependence of the differential cross section on the transverse momen- 
tum is more complicated. Indeed, for large pt polynomial decay is observed. 
For example, the ZEUS colloboration jl3| has fitted their measurements ob- 
tained in ep collision experiments by an empirical power law of the form 
(1 + const ■ Pt)' , with exponent a measured as a = 5.8 ± 0.5. 

In the following sections we will show that it is straightforward to ex- 
tend the Hagedorn theory in such a way that it decribes experiments at high 
energies (E > 10 GeV) as well. The basic tool for this is the recently de- 
veloped formalism of non-extensive statistical mechanics. It will just lead to 
asymptotic power laws of the above form. 

3 Tsallis' theory 

The formalism of non-extensive statistical mechanics is a generalization of 



the ordinary formalism of statistical mechanics ||14|-||16[. Wheras ordinary 



statistical mechanics is derived by extremizing the Boltzmann-Gibbs entropy 
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S = —J2iPi hipj (subject to constraints), in non-extensive statistical mechan- 
ics the more general Tsallis entropies 



*»~(i-I>?) 



(4) 



are extremized. The pi are probabilities associated with the microstates of 
a physical system, and q is the non-extensivity parameter. The ordinary 
Boltzmann-Gibbs entropy is obtained in the limit q — > 1. 

The Tsallis entropies are closely related to the Renyi information mea- 
|3T| and have similarly nice properties as the Boltzmann-Gibbs entropy 



sures 



has. They are positive, concave, take on their extremum for the uniform 
distribution, and preserve the Legendre transform structure of thermody- 
namics. However, they are not additive for independent subsystems (hence 
the name 'non-extensive' statistical mechanics). In the mean time a lot of 
physical applications have been reported for the formalism with q ^ 1. Exam- 



ples are 3-dimensional fully developed hydrodynamic turbulence |17|, [LJ| [L9 
2-dimensional turbulence in pure electron plasmas p3[ |24jl , Hamiltonian sys- 
tems with long-range interactions [SO], granular systems [ST], systems with 



strange non-chaotic attractors [p2] , and peculiar velocities in Sc galaxies |25| 
Given some set of probabilities Pi one can proceed to another set of prob- 
abilities Pi defined by 

' " (5) 



These distributions are called escort distributions [^ . Extremizing S q under 
the energy constraint 

£pe l = ^, (6) 

i 

where the e« are the energy levels of the microstates, the probabilities Pi 
come out of the extremization procedure as 

Pi = -^(1 + (9-1)^)" A, (7) 



where 



Z ? = £(1 + (g - l)/? ei )-^ (8) 



is the partition function and /3 — 1/T is a suitable inverse temperature vari- 
able (depending on U q ). In the limit q — > 1, ordinary statistical mechanics 
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is recovered, and the Pj just reduce to the ordinary canonical distributions 
Pj = ^e~ /3ei . For q ^ 1, on the other hand, they can be regarded as general- 
ized versions of the canonical ensemble, describing probability distributions 
in a complex non-extensive system at inverse temperature (3. 



4 Combining Hagedorn's and Tsallis' theory 
4.1 States of factorizing probabilities 

To generalize Hagedorn's theory we first have to decide on how to introduce 
grand canonical partition functions in non-extensive statistical mechanics. 
The problem is non-trivial, as one immediately recognizes from just consid- 
ering a system of N independent particles. In fact, it must first be defined 
what one means by independence in the non-extensive approach. The most 
plausible definition is that for independent particles probabilities should fac- 
torize. Let us consider such a distinguished factorized state where each joint 
probability is given by products of single-particle Tsallis distributions. This 
means, the joint probability to P^ia,...,^ to observe particle 1 in energy state 
e ix , particle 2 in energy state e i2 , and so on is given by 



1 N 



*W,i w = ^ IK 1 + (? - (9) 

At the same time, we could also describe our system by the Tsallis distribu- 
tion formed with the total energy (the Hamiltonian H(ii,t2, . . . , zV) of the 
system) 

/V„ , v = ^(1 + (<? - l)PH(h, • • • , i N T^- (10) 



Equating (g) and (|T0l) we see that the total energy of the system is given by 

N 

l + (q-l)(3H=l[(l + (q-l)(3e ij ), (11) 
i=i 

which we may write as 

H = E % + - + - E + ■■■ (12) 

3 3,k j,k,l 

(all summation indices are pairwise different). This means that for the unique 
particle state where probabilities factorize the total energy of the system is 
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not the sum of the single particle energies, provided q ^ 1. In other words, if 
seen from the energy point of view, the particles are formally interacting with 
a coupling constant (q — l)/3 although the probabilities factorize. This fact 
is not too surprising — the formalism of non-extensive statistical mechanics 
is of course designed to describe systems with (long-range) interactions, and 
also the entropy is non-additive. 

We can also invert the above statement. If we consider a Hamiltonian 
that is just the sum of the single particle energies, then the probabilities do 
not factorize provided q ^ 1. This is a well-known statement in non-extensive 
statistical mechanics. 

In the following, we will generalize Hagedorn's theory using the unique 
states of factorizing probabilities — thus implicitly introducing interactions 
between the particles from the energy point of view. 

4.2 Statistical mechanics of the fireball 

We will consider particles of different types and label the particle types by 
an index j. Each particle can be in a certain momentum state labelled by 
the index %. The energy associated with this state is 

€ij = Jpl + m? (13) 

(the relativistic energy- momentum relation). Let us define a non-extensive 
Boltzmann factor by 

Xij = (l + iq-l)^)'^. (14) 

It approaches the ordinary Boltzmann factor e _/3eiJ for q — > 1. We now very 
much follow Hagedorn's original ideas, replacing the ordinary Boltzmann 
factor by the generalized one. The generalized grand canonical partition 
function is introduced as 

y- I-:? as) 

O) ij 

Here denotes the number of particles of type j in momentum state %. The 
sum YLua stands for a summation over all possible particle numbers. For 
bosons one has = 0, 1, 2, . . . , oo, whereas for fermions one has = 0, 1. 
It follows that for bosons 

Y, x ^j = YZ: — (bosons) (16) 
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whereas for fermions 

^x^- 3 = 1 + Xij (fermions) (17) 



"%3 



Hence the partition function can be written as 

^ = nrr^Il( 1 + ^')> ( 18 ) 

ij L X V i'j' 

where the prime labels fermionic particles. For the logarithm of the partition 
function we obtain 

log Z = - £ Ml - Xij) + £ log(l + x vr ). (19) 

ij i'j' 

One may actually proceed to continuous variables by replacing 

H...H j£ -W [ ... 1 *-« J f rt . 1 * (20, 

{Vq\ volume of the interaction region) and 

£[•••] - / />H[-]^ (2i) 

where p(m) is the mass spectrum. 

Let us now calculate the average occupation number of a particle of 
species j in the momentum state i. We obtain 

^• = ^^-logZ = 



(l + (g-l)^)— ±1 



(22) 



where the — sign is for bosons and the + sign for fermions. 

In order to single out a particular particle of mass mo, one can formally 
work with the mass spectrum pirn) = 8(m — mo). To obain the probability 
to observe a particle of mass m in a certain momentum state, we have to 
multiply the average occupation number with the available volume in mo- 
mentum space. An infinitesimal volume in momentum space can be written 
as 

dp x dpydp z = (IplPt sin 9dpTd9 (23) 
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where = yp 2 + p 2 is the transverse monentum and p x = Pl is the longi- 
tutinal one. By integrating over all 9 and pi one finally arrives at a proba- 
bility density w(pt) of transverse momenta given by 

r°° 1 
w(p T ) = const ■ p T I dp L — . (24) 

J o (l + (q-l)djp 2 +p 2 +m 2 )—±l 



(1 + (q - l)P^p 2 T + p 2 L + m 2 ) «-i ± 1 
Since the Hagedorn temperature is rather small (of the order of the 7r mass), 



under normal circumstances one has j3 y p\ + p\ + m\ >> 1, and hence the 
±1 can be neglected if q is close to 1. One thus obtains for both fermions 
and bosons the statistics 

w(pr) ~ const ■ Pt J dp L (l + (q — l)f3\Jp^ + p\ + rnj)"^ 1 (25) 

which, if our theory is correct, should determine the px dependence of experi- 
mentally measured particle spectra. The differential cross section a~ l da /dpT 
is expected to be proportional to w(pt)- 



5 Large pt approximation 
5.1 The differential cross section 

To further proceed with analytic calculations, one may actually perform one 
further approximation step. This is the generalization of the step of going 
from eq. (|l]) to eq. ([3D in Hagedorn's original theory, which is a good approx- 
imation for large pt- Since for our applications q is close to 1, one expects 
that a similar step is possible in the more general non-extensive theory. 
Let us write the formula for the differential cross section in the form 

[ poo . , . 2_ 

- =cu dx (l + (q-l)Jx 2 + u 2 + m 2 g ) 9-1 . (26) 
a dpT Jo v v ' 

Here x = Pl/To, u = Pt/Tq and m@ := mo/T are the longitudinal mo- 
mentum, transverse momentum and mass in units of the Hagedorn tem- 
perature T , respectively, c is a suitable constant. Let us look at this 
formula for large values of pr- If u is very large, we can approximate 
^Jx 2 + v? + m 2 p = uyl + (x 2 + m 2 ^)ju 2 ~ u + (x 2 + m 2 ^)/(2u). Of course, 
for this to be true the integration variable x should not be too large, but for 
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large x the integrand is small anyway and yields only a very small contribu- 
tion to the cross section. Moreover, since u is large we may also neglect the 
mass term m 2 g, arriving at the aprroximation 




<73) 

— q 

1 drr r°° I I r 2 \ \ 9-1 

/ . / . / .1 \\ (27) 

Although this may look like a rather crude approximation, in practice it is 
quite a good one. This is illustrated in Fig. 1, which shows the right hand 
sides of eq. (p6|) and eq. ( p?D for q = 1, 1.1, 1.2. The lines corresponding to 
the exact expression (|26|) and the approximation (|27D can hardly be distin- 
guished. The range of u and q is similar to what we will use later for the 
comparison with experimental measurements. Since the Hagedorn tempera- 
ture T is about 120 MeV, and since very often pions are produced, we have 
chosen in formula fl26|) m 2 g = (m^/To) 2 ~ 1.3. The errors made in eq. (|27|) 
by neglecting mg and by approximating the square root work in opposite 
directions and almost cancel each other for q « 1.1. 

Within this approximation we can now easily proceed by analytical means. 



We may write eq. (|27|) as 



C. i_ 
1 + - — j- — , —^x 2 28 
2u(l + (q - l)u) J v ' 

Substituting 

' !/ ~ 1 x (29) 



V 2u(l + (q- l)u) 
and using the general formula p3| 

roo t 2x ~ l 1 



poo t | 

I (i + t 2 y+y dt = 2 B{ ~ x ' y) ( Rex >^ Re y>°)> (3°) 



where 



T(x + y) 

denotes the beta-function, one arrives at 



B (,, y ) = ™> (31) 



1 da 
a dpT 



c(2(q - 1))-V 2 5 (I, - ^ u 3/ 2 {1 + {g _ l)u y^ + h . 



(32) 



This formula, with suitably determined c, q, T , will turn out to be in very 
good agreement with experimentally measured cross sections. 



10 



5.2 Normalization 

Let us quite generally consider a probability density with the above it-depen- 
dence 

p(u) = — u 3/2 (l + (q- (33) 
Normalization yields 

Z* 00 q 1 

Z q = u 3/2 (1 + (q - l)uy— + z du. (34) 

J 



Substituting t 2 := (q — l)u the integral can be brought into the form ( |30l) 
and one obtains 



Zi = (9 _ 1) -v. B g,_±__,) 



(35) 



Note that the beta function further simplifies if — 3 is an integer. For 
n £ N one generally has 

(n - IV 

5(x,n) = - w 1 . J \ r, (36) 



in our case x = |. 



5.3 Moments 

Let us now evaluate the moments of w defined by 

/•oo 1 /*oo o a 1 

= / u m p(u)du = — / n5+ m (l + (q- l)u)-— + * du. (37) 
jo Z q Jo 

Again substituting t 2 = (q — l)u the integral can be evaluated to give 

— 3 — ml 

(38) 

Generally the Gamma function satisfies 

T(x + l)=xT(x) (39) 

11 



i r (| + m) 


r 




— mj 


( q -i) m r(l) 


r 







which one may iterate to obtain 



m—1 

T(x + m) = T(x) J| (x + j). (40) 

3=0 



Using this in eq. (|38|) one finally arrives at 



1 m_1 ^ 4- 9i 

t u ™\ = — TT J . (41) 



In particular, one obtains for the average of u 



2 4 -3q 

and for the second and third moment 



(u 2 ) = -— — (43) 

W 44-3g5-4g v ; 

^ ~ 84-3g5-4g6-5g' ^ 



The variance is given by 



" 2 ■= <""> - ^ ! (4 -3^(5-4,) - < 45 » 
Note that the moments obey the simple recurrence relation 

(u m+1 ) = (u m ) ■ A4 _ 5+ f™ ■ (46) 
4 + m — (3 + m)q 

6 (/-dependence of the Hagedorn temperature 

A fundamental property of Hagedorn's theory is the fact that the Hagedorn 
temperature T is independent of the center of mass energy E. Now, in the 
generalized theory we have a new parameter q and it is a priori not clear if 
and how T depends on q. 

However, looking at eq. ( |T2"D we recognize that if q increases from 1 to 
slightly larger values the (formal) interaction energy of our system increases. 
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This increase in interaction energy must come from somewhere and is ex- 
pected to be taken from the (finite volume) heat bath of the fireball. Thus 
the effective temperature of the bath is expected to slightly decrease with 
increasing q. 

To get a rough estimate of this effect, let us work within the approxima- 
tion scheme of the previous section. Differentiating eq. ([33|) with respect to 
u one immediately sees that the distribition p{u) has a maximum at 

"* = | = 3^' < 47 > 

On the other hand, the experimentally measured cross sections always ap- 
pear to have their maximum at roughly the same value of the transverse 
momentum p?, namely at 

p* T « 180 MeV, (48) 

independent of the beam energy E. This implies that in the generalized 
thermodynamic approach the effective Hagedorn temperature To will become 
(slightly) g-dependent. 

To = (l - |) P* T (49) 

The variation of T with energy E is actually very weak. In the next section 
we will consider two extreme cases, namely q — > 1 for E — > and q — > ^ for 
E — > oo. The first case corresponds to T = 120 MeV, the second case to 
To = 107 MeV. This is only a very small variation. 

7 Energy dependence of q 

We still have to decide how the non-extensitivity parameter q depends on 
the center of mass energy E of the beam. We will present some theoretical 
arguments, which are well supported by the experimental data. 

7.1 Plausible value for q max 

Let us define a to be the power of the term (1 + (q — 1)«) _1 in eq. (|33|) 

a = i=i - \ < 50) 
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Clearly, for small energies E Hagedorn's theory (q = 1) is valid: 

E -> a -> oo g -> 1 (51) 

For increasing E 1 , a should decrease. For example, the ZEUS collaboration 
measures a = 5.8 ± 0.5 at medium energies ||13|| . However, a cannot become 
arbitrarily small, because we must have a finite variance of u (otherwise 
statistical mechanics does not make sense). For large u, we can certainly 
use the large pr approximation of section 5. Asymptotically the density p(u) 
decays as 

p{u) ~ u - a+ ^ 2 (52) 
Thus u 2 p(u) decays as u~ a+7 l 2 and hence (u 2 ) = J du u 2 p(u) only exists if 

7 9 5 

-a + -<-l^a>-^q<- (53) 

(for a = 9/2 there is a logarithmic divergence). This is a strict upper bound 
on q. 

If, in addition, we postulate that the distinguished limit theory corre- 
sponding to the largest possible value of q should be analytic in y/u, then 
only integer and half-integer values of a are allowed. This leads to 

aw = 5 q max = — = 1.222 (54) 

y 

as the smallest possible value of a, respectively largest value of q. 



7.2 A smooth dependence on q 

It is reasonable to assume (and supported by the experimental results) that 
the smallest possible value of the exponent a is taken on for the largest 
possible energy E. Hence we have the two limit cases 

E -> a -> oo (55) 

E -> oo =^ a -> a min = 5. (56) 

A smooth monotonously decreasing interpolation between these limit cases 
is given by the formula 
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where E is some suitable energy scale. It turns out that the experimental 
data are perfectly fitted if we choose 



E = -m z = 45.6 GeV, (58) 

mz being the mass of the Z° boson. 

Whether this coincidence of the relaxation parameter E with half of the 
Z° mass is a random effect or whether there is physical meaning behind this 
is still an open question. What is clear is the following. At center of mass 
energies of the order 90 GeV often a massive elektroweak gauge boson Z° is 
produced. If a Z° at rest decays into a quark pair qq, then each quark has 
the initial energy \mz, which is input to the hadronization cascade. It is 
obvious that at this energy scale there is a rapid change in the behaviour of 
the fireball due to the occurence of the electroweak bosons. The same mass 
scale now seems to mark the crossover scale between ordinary Hagedorn 
thermostatistics, valid for E « \mz, and generalized thermostatistics with 
q ~ qmax, valid for E » \m z . 

Eq. (]57|) can equivalently be written as 

11 _ e -E/E 

q(E) = (59) 



8 Energy dependence of the multiplicity 

The differential cross section a~ l daj dpx measured in experiments is not a 
normalized probability density. This is due to the fact that it is dependent 
on the multiplicity M (the average number of produced charged particles). 
Also, it has dimension GeV" 1 , whereas the probability density p(u) is di- 
mensionless. All this, however, is just a question of normalization. 

To proceed from p(u) to a~ x da /dp?, we should multiply p{u) with the 
multiplicity M. Also, in order to give the correct dimension of a cross section, 
we should multiply with T _1 . Thus we arrive at the formula 

1 da 1 , , „ . „ 

-— = -Mp (u . 60 
a dp T T 

The multiplicity M as a function of the beam energy E has been indepen- 
dently measured in many experiments. A good fit of the experimental data 
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in the relevant energy region is the formula 

/ E \5/H 

M = — T To 9 " 1 = 120 MeV (61) 



see Fig. 2, data as collected in [34|). It is remarkable that the constant in 



front of this power law is 1 if E is measured in units of the Hagedorn Tem- 
perature Tq -1 . Apparently, the Hagedorn temperature is a very appropriate 
portion of energy for our statistical approach. Essentially, the scaling law 
([61] ) says that only a certain fraction of the logarithm of the energy E (de- 
scribed by the scaling exponent 5/11 ~ 0.45) is used to increase the number 
of charged particles. 



9 Comparison with experimentally measured 
cross sections 

We are now in a position to directly compare with experimental measure- 
ments of cross sections. All parameters of the generalized Hagedorn theory 
(the Hagedorn temperature T (q), the non-extensitivity parameter q(E) and 
the normalization constant of the cross section) have been discussed in the 
previous sections and concrete equations have been derived. 

Formula ( |60D with p{u) given by eq.(^3|), q(E) given by eq. (|59l), T (q) 
given by eq.(|4^) and multiplicity M(E) given by eq. (|6"T| ) turns out to very 
well reproduce the experimental results of cross sections for all energies E. 
This is illustrated in Fig. 3, which shows experimental cross sections mea- 
sured by the TASSO and DELPHI collaborations (|Tl], [L^], data as collected 
in ||34|| ) as well as our theoretical prediction. For all curves we have chosen 
the same universal parameters p* T = 180 MeV and E = \rnz-, so we do not 
use any energy- dependent fitting parameters. 

The agreement is remarkably good. Hence the statistical approach pre- 
sented in this paper qualifies as a simple thermodynamic model that explains 
the experimental data quite well. In particular, for the first time analytical 
formulas are obtained that correctly describe the measured cross sections. 
The quality of agreement is at least as good as that of Monte Carlo simula- 
tions, in spite of the fact that Monte Carlo simulations usually use a large 
number of free parameters. For small pt, the agreement seems even to be 
slightly better than that of Monte Carlo simulations (see e.g. |34j for typical 
Monte Carlo results). 
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Our approach yields g-values of similar order of magnitude as the ones 
used in the fits in ||29|| . In our theoretical approach all parameters are 
now given by concrete formulas. Some of these formulas, most importantly 
eq. (0), have been derived empirically and should thus be further checked 
and possibly refined by further experimental measurements. 

10 Conclusion 

In this paper we have developed a thermodynamic model describing the 
statistics of transverse momenta of particles produced in high-energy col- 
lisions. Our approach is based on a generalization of Hagedorn's theory 
using Tsallis' formalism of non-extensive statistical mechanics. Hagedorn's 
original theory is recovered for small center of mass energies, whereas for 
larger energies deviations from ordinary Boltzmann-Gibbs statistics become 
relevant. The crossover energy scale between the two regimes E — > 0, q ~ 1 
and E — > oo, q m ^ is given by half of the Z° mass. 

At large energies the generalized thermodynamic theory implies that cross 
sections decay with a power law. This power law is indeed observed in exper- 
iments. Our theory allows us to analytically evaluate moments of arbitrary 
order (within a certain approximation). We obtain formulas for differential 
cross sections that are in very good agreement with experimentally measured 
cross sections in annihilation experiments. 

We suggest that future analysis of the experimental data should concen- 
trate on precision measurements of the energy dependence of the parameter 
q. For this the mass spectrum of produced particles should carefully be taken 
into account. Generally, it appears that high-energy collider experiments do 
not only yield valuable information on particle physics, but they may also be 
regarded as ideal test grounds to verify new ideas from statistical mechanics. 
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Figure captions 



Fig. 1 Comparison between the exact formula ( p6[ ) with m/| = (m^/To) 2 = 1.3 
and the approximation (|27| ) for various values of q. The constant c was set equal 
to 1 for all q. 

Fig. 2 Dependence of the multiplicity M on the center of mass energy E. The 
figure shows the experimental data and a straight line that corresponds to the 
scaling law (El]). 

Fig. 3 Differential cross section as a function of the transverse momentum pt 
for various center of mass energies E. The data correspond to measurements of 
the TASSO (E < 44 GeV) and DELPHI [E > 91 GeV) collaboration. The solid 
lines are given by the analytic formula (|60|). 
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